Self-localization of a small number of Bose particles in a superfluid Fermi system 
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We consider self-localization of a small number of Bose particles immersed in a large homogeneous 
superfluid mixture of fermions in three and one dimensional spaces. Bosons distort the density of 
surrounding fermions and create a potential well where they can form a bound state analogous to a 
small polaron state. In the three dimensional volume we observe the self-localization for repulsive 
interactions between bosons and fermions. In the one dimensional case bosons self-localize as well 
as for attractive interactions forming, together with a pair of fermions at the bottom of the Fermi 
sea, a vector soliton. We analyze also thermal effects and show that small non-zero temperature 
affects the pairing function of the Fermi-subsystem and has little influence on the self-localization 
phenomena. 

PACS numbers: 03.75.Ss, 03.75.Hh, 03.75.Lm 
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I. INTRODUCTION 

Ultra-cold atomic gases offer possibilities for realiza- 
tions of complex mathematical models used in different 
fields of physics with an unprecedented level of the ex- 
perimental control For example, condensed matter 
phenomena like the superfluid-Mott insulator transition 
and the Bose-glass phase or the Anderson localization ef- 
fects can be experimentally investigated . Fermionic 
gases, in particular Fermi superfluids, have received a lot 
of attention, especially after the observation of the tran- 
sition between the superfluid Bardeen-Coopcr-Schrieffcr 
(BCS) pairs and the Bose-Einstein condensate (BEC) of 
diatomic molecules 0, Q . 

The behavior of a small object immersed in degenerate 
quantum gases has been investigated by several authors 
[9l-[l8|. For example, weak interactions between a sin- 
gle impurity atom and particles of a large BEC can be 
described by the perturbation theory. For stronger inter- 
actions an effective mass of an impurity atom diverges 
indicating the breakdown of the perturbation approach 
and the self-localization of the impurity object in a close 
analogy to the small polaron problem, i.e. localization of 
an electron in a surrounding cloud of lattice distortions 
[igj . In ultra-cold fermionic gases an example of polaron 
effects with a small number of spin- up fermions immersed 
in a large cloud of spin-down Fermi particles has been 
studied theoretically j2(]|425| and recently realized exper- 
imentally [2^[23|. Employing a Feshbach resonance, that 
allows tuning the interaction strength between atoms, ex- 
perimentalists have been able to investigate a transition 
from the nearly non-interacting case, through the polaron 
regime to the limit where pairs of unlike fermions form 
tightly bound molecules. 

In the present publication we consider a small number 
of Bose particles immersed in a large, homogeneous, su- 
perfluid and balanced mixture of spin-up and spin-down 
fermions and analyze the self-localization phenomenon. 
Another limit, investigated already in the literature, con- 
cerns Bose- Fermi mixtures with a number of bosons com- 



parable to (or even larger than) a number of fermions and 
effects of the phase separation [28l - l38j . The latter cor- 
responds to instability of a homogeneous solution when 
boson-fermion interaction reaches a critical strength. In 
the case of small boson numbers, the boson-boson inter- 
actions can be neglected and the uniform density solution 
is unstable as soon as the boson-fermion coupling con- 
stant becomes non-zero. However, this does not mean the 
self-localization of Bose particles. We show that the self- 
localization takes place for stronger interactions when the 
boson-fermion coupling constant is greater than a non- 
zero critical value. 

The possibility of solitonic behavior in Bose- Fermi mix- 
tures with fermions both in the normal and superfluid 
states has been investigated in the literature [33, 39 41j]. 
For a large number of bosons, if the attractive boson- 
fermion interaction is sufficiently strong, the boson-boson 
repulsion may be outweighed and the whole Bose and 
Fermi clouds reveal solitonic behavior. We consider Bose- 
Fermi mixtures in the opposite limit of small boson num- 
bers. In that regime different kind of solitons exists. In- 
deed, in the ID case description of the system may be 
reduced to a simple model where bosons and a single pair 
of fermions at the bottom of the Fermi sea are described 
by a vector soliton solution. 

The paper is organized as follows. In Sec. |ll] we in- 
troduce the model used in the description of Bose- Fermi 
mixtures. The results for the case of three-dimensional 
(3D) and ID spaces are collected in Sec. Illll and we con- 
clude in Sec. IIVI 



II. MODEL DESCRIPTION 

Let us consider a small number TVf, of bosonic atoms in 
the Bose-Einstein condensate state immersed in a homo- 
geneous, dilute and balanced mixture of fermions in two 
different internal spin states in a 3D volume. Interactions 
of ultra-cold atoms can be described via contact poten- 
tials Vij{r) = gij6{r) with strengths given in terms of 
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where m, 



s-wave scattering lengths aij as gij 
stands for a reduce mass of a pair of interacting atoms. 
In our model we consider attractive interactions between 
fermions in different spin states, i.e. negative coupling 
constant gff. Interactions between bosons and fermions 
are determined by the spin- independent parameter g^f. 
We neglect mutual interactions of bosonic atoms in the 
assumption that either their density remains sufficiently 
small or the coupling constant is negligible. 

The system is described by the following Hamiltonian 



H = 



dV 
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where i/o = 



^ — /i. ^'b and ^ f^s refer, respectively, 
to the field operators of bosonic and fermionic atoms 
where s € {+,— } indicates a spin state, fi stands for 
the chemical potential of the Fermi sub-system and rrib 
and ruf are masses of bosons and fermions, respectively. 

We look for a thermal equilibrium state assuming that 
the Bose and Fermi sub-systems are separable. For in- 
stance in the limit of zero temperature it is given by a 
product ground state 
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We also postulate that the Fermi sub-system can be 
described by the BCS mean-field approximation Q 

with the paring field A(r) = \gff\ i'lpf^+'ilif^-j and the 



Hartree-Fock potential W{y) 



9ff\ \^f,+^f 



\9ff \ -''Pf,-^ affected by a potential proportional 
to the density of bosons A^f,|(/)(r)p. Assuming a spheri- 
cal symmetry of particle densities, the description of the 
system reduces to the Bogoliubov-de Gennes equations 
for fermions 

+ ^Vnlm = Enl 
^*Unlm - {Hq + W + gbfNb\(j)\'^) Vnlm = EnlVnlm, 

(3) 

where / and m stand for angular momentum quantum 
numbers and 

W = -|5//lE[^"'l^"'™WI' + (l-^"')l^"'™WI']' 

nlm 

(4) 

A = |g//|E(l-2/«0«»;™(r)<z™(r), (5) 
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with the Fermi-Dirac distribution 
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which have to be solved together with the Gross- 
Pitaevskii equation for bosons 
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W{r) = gbfPfir). 
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The effective potential V{r) for bosons comes from con- 
tact interactions between bosons and fermions. P/ is 
density of fermions and fib is the chemical potential for 
bosons. We consider the temperature much lower than 
the critical temperature for Bose-Einstein condensation 
therefore we can neglect thermal excitations of bosons. 

The coupled equations ([3]) and (O are solved numer- 
ically in a self-consisted manner. In the calculations we 
adopt 
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units for energy and length, respectively, where kp = 
(STr^no)^/^ is the Fer mi wave-number of a uniform ideal 
Fermi gas of density uq. In these units the coupling con- 
stants are the following 
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and we deal with six independent parameters in the sys- 
tem: number of bosons TV;,, chemical potential of Fermi 
sub-system ^, ratio of the masses scattering lengths 
kpaff and kpabf and radius R of the 3D volume we con- 
sider. 

In the 3D case the coupling constant gff in A [Eq. ([5])] 
has to be regularized in order to avoid ultraviolet diver- 
gences. That is, gff — > g^s where 
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The logarithmic term in (jlip results from the sum over 
Bogoliubov modes corresponding to the energy above a 
numerical cutoff Ec performed in the spirit of the local 
density approximation, see j42| - [45| for details. 



III. RESULTS 

Without interactions between bosons and fermions the 
ground state of the system corresponds to uniform par- 
ticle densities. For the non-zero coupling constant gbf, 
the uniform solution become unstable and, depending on 



Figure 1: (Color online) Self-localization of '^^Na atoms in a 
superfluid mixture of ''"K atoms. Panel (a) shows the pair- 
ing function A(r), panel (b) fermion density pf{r) and panel 
(c) density of bosons \(j){r)\'^. Solid black lines correspond 
to boson- fermion interaction strength gbf = 10 and dotted- 
dashed blue lines to gtf — 0. In panel (c) the dotted-dashed 
blue line is hardly visible because for gi,f — bosons are de- 
localized and their density very small. Number of bosons 
Nt = 100 and fermions Nf ~ 12000 (chemical potential 
fi — Ef) and fermion- fermion coupling constant gff — —5.5. 



the sign of gi,f, the bosonic and fermionic clouds tend 
to separate from each other or try to stick together. 
For sufficiently strong interactions, the effect of the self- 
localization may be expected (see the similar problem in 
the case of an impurity atom immersed in a large Bose- 
Einstein condensate considered in Ref. [9l-[T]|). Indeed, 
for gtf > bosons repel fermions and create a potential 
well in their vicinity where they may localize if the well 
is sufficiently large. For attractive interactions the den- 
sity of fermions increases in the vicinity of Bose particles. 
Due to the fact that gi,f < 0, the bosons experience the 
density deformation in a form of a potential well and they 
may localize. 

We begin with the 3D model and focus on the repul- 
sive boson-fermion interactions. Analysis of both zero- 
temperature limit and thermal effects are performed. 
Then we consider the ID case where the self-localization 
phenomenon may be related to the presence of a vector 
soliton solution. 



Three dimensional model 



Figure 2: (Color online) Probability densities |Wnim('')P of 
two fermion pairs at the bottom of the Fermi sea with angular 
momentum I = 0. Panel (a) corresponds to the ground state 
(n = 1) of the radial degree of freedom and panel (b) to 
the first excited state (n = 2). Solid black lines correspond 
to boson-fermion interaction strength gbf — 10 and dotted- 
dashed blue lines to gbf — 0. All parameters are the same as 
in Fig. m 
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Figure 3: (Color online) The average radius of the Bose cloud 
(r) [panel (a)] and the standard deviation a = {{r^} — {r)^)^^^ 
[panel (b)] versus boson-fermion coupling constant gbf. Black 
full symbols correspond to a mixture of ^^Na and '''^K atoms 
while red open symbols to a mixture of ^Li and ^''K atoms. 
Note the abrupt transitions to localized states when critical 
values of gbf are reached. All the other parameters are the 
same as in Figs. [TJ21 



Figure [T] shows the densities of bosons and fermions 
and the pairing function corresponding to the ground 
state of the system for gi,f = and gbf = 10. With- 
out boson-fermion interactions the quantities are flat and 
uniform (except a small region close to the edge of the 
3D volume due to assumed open boundary conditions). 
However, when the considerable interactions are turned 
on it becomes energetically favorable to separate bosons 
and fermions, the p/(r) is depleted around the center and 
bosons form a bound state localized in small area around 



r = 0. It is clear, that the localization effect is the result 
of boson-fermion interactions. It relies on a local defor- 
mation of the density of fermions and is not affected by 
the boundary conditions. 

The response of the Fermi sub-system to bosons, that 
tend to localize, can be investigated by monitoring de- 
formation of the Bogoliubov quasi-particle modes. The 
density of fermions is the sum of the Bogoliubov modes 
Pfi^) ='^J2nim l^'n/m(r)P- The modcs with zero angular 
momentum contribute only to the density around r = 0. 
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Figure 4: (Color online) Self-localization of '^^Na atoms in a 
superfluid mixture of '*''K atoms for non-zero temperature. 
Panel (a) shows the pairing function A(r), panel (b) fermion 
density p/(r) and panel (c) density of bosons |0(r)p. Solid 
black lines correspond to T = and ^ = Ef, red dashed lines 
to T = 0.028Tf and fi = Ef, blue dotted-dashed lines to T = 
6Tf and fi = 0.16Ef- Boson-fermion interaction strength 
Qbf ~ 10, fermion-fermion coupling constant gff — —5.5 and 
number of bosons A^;, = 100 and fermions Nf ~ 12000. In 
panel (a) the dotted-dashed blue line is not visible because 
for T = 6Tf the pairing function is equal zero. In panels 
(c) and (d) the solid black and dashed red lines are hardly 
distinguishable . 



Consequently, the modification of these modes is primar- 
ily responsible for preparation of the potential well in 
which bosons localize. In Fig.[3]we illustrate the deforma- 
tion of two modes with / = corresponding to fermions 
at the bottom of the Fermi sea but wc should keep in 
mind that all modes with / = become affected by the 
interactions with bosons. The deformation of modes for 
fermions at the Fermi level is reflected by a change of a 
shape of the pairing field visible in Fig. [TJ because those 
modes contribute mainly to A(r). 

The interaction of fermions and the impurity Bose par- 
ticles influences the pairing function A only locally, see 
Fig. m It implies that the superfluidity is not destroyed 
even when the interaction is so strong that the localiza- 
tion of the impurity object takes place. 

The data in Figs. [T][l] are related to Nb = 100 ^^Na 
atoms and the mixture of Nf « 12000 ^'^K atoms (chem- 
ical potential /i = Ep) in two different hyperfine states. 
We set the scattering lengths gff = —5.5 and gtf = 10 
with the assumption that they can be realized by the 
use of the Feshbach resonances (e.g. magnetic resonance 
for fermions and optical resonance between bosons and 
fermions H^). In Fig. [3] we show the average ra- 
dius of the Bose cloud (r) and the standard deviation 
a — {{r^) — {r)^y/^ as a function of the coupling con- 
stant gbf. The self- localization means that both (r) and 
a are much smaller than the radius of the 3D volume. 
One can see that there is a critical non-zero value of gtf 
leading to the self- localization. This critical gi,f is dif- 
ferent from the critical value for the instability of the 



homogeneous solution (i.e. phase separation condition). 
The latter, for the case without boson-boson interactions, 
corresponds to gt/ > 0. If we replace the sodium atoms 
by ^Li atoms, it turns out that the critical value of gbf 
for the self-localization increases. This is, because com- 
pressing the cloud of light lithium particles costs more 
energy than in the case of heavier sodium atoms. 

A small non-zero temperature mostly affects superflu- 
idity and has little effect on the self-localization phe- 
nomenon. Indeed, in Fig. U] we see that even for T = 
0.028rF when the pairing function is very small the den- 
sities of bosons and fermions hardly change. Increas- 
ing temperature to T = Tp (which is still much smaller 
than critical temperature for Bose-Einstein condensation 
of Nb ~ 100 bosonic atoms localized in a volume of the 
radius (r) w 4, i.e. Tbec ~ GTp) wc observe effects of 
thermal fluctuations in the fermion density and a modi- 
fication of the density of bosons but the self-localization 
persists. Thus, bosons self- localize both for the normal 
and superfluid phase of the Fermi sub-system. 

We have considered the repulsive boson-fermion inter- 
action. For the attractive interaction we do not observe 
the self-localization regardless on the phase of the Fermi 
sub-system. For gbf < the particle densities may col- 
lapse to Dirac-delta distributions. For sufficiently small 
\gbf\ a metastable state may appear. However, it turns 
out that the existence of such a metastable state is not 
the result of self-localization in the system. Indeed, it 
is an effect of a compromise between the requirement of 
minimal kinetic energies and restrictions related to the 
boundary conditions. In the following we consider a ID 
model where there is no problem with the collapse of the 
densities and show that Bose particles can localize in the 
Fermi sub-system for attractive boson-fermion interac- 
tions too. 



B. One dimensional model 

If in X and y directions we apply harmonic potentials 
of frequency ui± and huj± is much greater than the chem- 
ical potentials, the 3D system becomes effectively one- 
dimensional. Then, the description reduces to the ID 
version of Eqs. ©-([HI) with the following coupling con- 
stants [in the units ©] 
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(12) 



which have been obtained assuming that the x and y de- 
grees of freedom of each atom are in the ground states of 
the harmonic potentials. In the ID case there is no ul- 
traviolet divergence and the pairing function does not re- 
quire regularization. Nevertheless numerical simulations 
converge much faster if the Bogoliubov modes, above a 
numerical cut-off energy Ec, are included in the spirit of 
the local density approximation. That is, the coupling 
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Figure 5: (Color online) Self-localization of a single boson 
(A^i, = 1) in a superfluid mixture of fermions in ID space. 
Panel (a) shows the pairing function A(z), panel (b) fermion 
density Pf(z) and panel (c) boson density |(^(z)|^. Solid black 
lines correspond to boson- fermion interaction strength gtf ~ 
—20 and dotted-dashed blue lines to Qbf ~ 0. Number of 
fermions Nf ~ 20 (chemical potential fi = Ef) and fermion- 
fermion coupling constant gff — —1. Ratio of masses of Bose 
and Fermi particles fulfills Eq. (|16|l . The configuration 

space extends from z = — 10 to 2 = 10. In panel (c) the 
dotted-dashed blue line is hardly visible, because the boson is 
delocalized and its density very small for gtf = 0. Red dashed 
line in panel (c) indicates the solitonic solution Eq. p7p . 
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Figure 6: (Color online) Bogoliubov modes Wfe(z) correspond- 
ing to fermion pairs located close to the bottom of the Fermi 
sea. Panel (a) is related to the pair of fermions at the bottom 
of the Fermi sea, panel (b) and (c) to the next pairs. Solid 
black lines correspond to the numerical solutions. Red dashed 
line in panel (a) indicates solitonic solution Eq. (|17p . All the 
others parameters are the same as in Fig. [S] 



constant in A is substituted by 
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For repulsive boson-fermion interactions, we observe 
the self-localization of bosons with the behaviour of the 
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Figure 7: (Color online) Width of the boson density, i.e. 
(7 = ((z^) — (2;)^)^''^, versus boson-fermion coupling constant 
Qiyf. Black full symbols correspond to the numerical values 
and red open symbols to the solutions Eq. (|17|) . The configu- 
ration space extends from z = — 20 to z = 20. All the others 
parameters are the same as in Fig. (5] 



particle densities similar as in the 3D case. Therefore we 
focus on attractive interactions only. Figure [5] shows the 
results for g?,/ = ^20, obtained with periodic boundary 
conditions for fermions and open boundary conditions 
for bosons. For the attractive interactions bosons and 
fermions try to stick together which leads to an increase 
of the fermion density in the vicinity of the boson concen- 
tration and the creation of a potential well for localization 
of Bose particles. 

Analyzing the BogoHubov modes Vk{z) (see Fig. |6]) 
we find out that the probability density v^^z) of a pair 
of fermions at the bottom of the Fermi sea becomes 
strongly localized. The Bogoliubov mode vi{z) of the 
next fermion pair forms also a bound state. Since Vi{z) 
is an antisymmetric function it is nearly zero in the area 
around z = 0. Probability densities of other fermions 
arc deformed and almost all of them drop to zero in the 
region where ^0(2;) is localized. This may be interpreted 
as a realization of the Pauli exclusion rule. In the BCS 
limit only particles close to the Fermi level contribute to 
the pairing function A and there is practically no contri- 
butions from fermions located deeply in the Fermi sea. 
Therefore there is also no contribution from the pair of 
fermions at the bottom of the Fermi sea. That is why 
A(z), contrary to the fermion density, reveals a minimum 
at z = 0, see Fig. [SJ 

The analysis of the Bogoliubov modes suggests a sim- 
ple model of self-localization in the case of attractive 
boson-fermion interactions. Suppose, that in the vicinity 
of the localized bosons we may neglect the pairing field 
and the density of all fermions except a fermion pair at 
the bottom of the Fermi sea. Then, we obtain the fol- 
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lowing set of equations 
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there exists analytical solution of Eqs. p4|)-(fT5 
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Such a solution resembles vector solitons. They ap- 
pear in non-linear optics when interactions of several 
field components are described by a set of coupled non- 
linear Schrodinger equations 47[. Note that for the self- 
localization of an impurity atom in a large BEC consid- 
ered in Ref. [ll|, the ID system is described by a para- 
metric soliton with the state of the impurity atom given 
by the hyperbolic secant squared function. 

A comparison of the analytical solutions (|17p with nu- 
merical results of the full set of equations is shown in 
Figs. [5][6l The agreement is very good and increases 
with the strength of boson-fermion interactions. Indeed, 
for the strong interaction, due to the Pauli exclusion 
rule, there is negligible probability density to find other 
fermions than the localized pair in the vicinity of z = 0. 
As a consequence, the localized bosons interact almost 
exclusively with the localized fermion pair and the set of 
Eqs. (fT4l) -(fT5 |) becomes exact. 

Figure[6lD shows that the Bogoliubov mode vi (z) forms 
an antisymmetric bound state. In the vicinity of z = 
(where the fermion density is dominated by Vq and the 
pairing function drops to zero) this mode should fulfill 
equation similar to Eq. ([T4|) . that is 
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If (j) and I'o are given by Eq. (jl7p the antisymmetric so- 
lution of Eq. (|19p forms a marginal bound state 
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In the full description of the system, the state govern 
by the equation (19) may become cither truely bound or 
unbound. In the considered system, it turns out that the 
state is pushed towards a true bound state as visible in 
Fig. 6b. 

When the boson-fermion coupling constant gbf is de- 
creased, we observe the increasing discrepancy between 
analytical and numerical solutions, see Fig. [71 The width 
of the boson probability density obtained numerically 
is significantly greater than the corresponding analyti- 
cal value. This is due to the fact, that in the effective 
potential experienced by the bosons a considerable con- 
tribution comes from other fermions, and not only from 
the pair at the bottom of the Fermi sea. The density of 
such fermions, contrary to the localized fermion pair, pos- 
sesses a minimum at z = and thus effectively makes the 
potential well for bosons weaker. Consequently, bosons 
occupy a much larger space than can be expected on the 
basis of solutions Eq. ([T7| . 



IV. CONCLUSIONS 

We have considered a small number of bosons im- 
mersed in a superfiuid mixture of fermions in two differ- 
ent spin states. With negligible boson-boson interactions, 
homogeneous densities of the particles become unstable 
as soon as the boson-fermion coupling constant is non- 
zero. It corresponds to the phase separation transition. 
We show that in 3D space for sufficiently strong repul- 
sive boson-fermion interactions another transition takes 
place, i.e. the self-localization of Bose particles. That is, 
the repulsion between particles creates a local potential 
well for bosons where, if the well is sufficiently large, they 
can localize. The self-localization is present both for the 
superfiuid and the normal state of fermions. It modi- 
fies properties of the Fermi sub-system locally without 
destroying the superfiuidity. Low non-zero temperature 
affects the pairing function but has little effect on the 
self-localization phenomenon. 

We do not observe the self-localization for attractive 
boson-fermion interactions in the 3D case. In this context 
the self-localization requires sufficiently strong boson- 
fermion interactions. However, for strong attractive in- 
teractions no metastable state of the system exists and 
the densities of the atoms collapse to Dirac-delta distri- 
butions indicating a breakdown of the description with 
the contact interaction potentials. In the ID case there 
is no collapse for attractive boson-fermion interactions. 
The self-localization of bosons is accompanied by local- 
ization of a pair of fermions at the bottom of the Fermi 
sea. This phenomenon can be described by a simple 
model where the self-localization is related to the exis- 
tence of a vector soliton solution. 

To realize experimentally the self- localization of bosons 
in a Fermi system, ultra-cold clouds of bosons and 
fermions have to be prepared in a laboratory. For a 
sufficiently large boson-fermion coupling constant, that 
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can be achieved by means of a Feshbach resonance, 
the self-locahzation takes place. Signatures of the self- 
locahzation can be visible in expansion of the atomic 
clouds after trapping potentials are turned off. That 
is, if during the time of flight the boson-fermion inter- 
actions are kept negligibly weak, the initially strongly 
localized boson cloud will show much faster expansion 
than the Fermi cloud due to release of a large kinetic 
energy. The simplest experiment would employ a Fermi 
sub-system in a normal phase. In order to observe the 
self-localization in a superfluid Fermi mixture a manip- 
ulation of a fermion-fermion coupling constant is also 



needed and two Feshbach resonances must be employed, 
e.g. one resonance controlled by magnetic field and the 
other by optical means. 



Acknowledgments 

This work is supported by the Polish Government 
within research projects 2009-2012 (KT) and 2008-2011 
(KS). 



M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, [22 
A. Sen, and U. Sen, Adv. Phys. 56, 243 (2007). 
I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. [23 
80, 885 (2008). 

M. Greiner, O. Mandel, T. Esslinger, Th. W. Hansch, [24 
and I. Bloch, Nature 415, 39 (2002). 

L. Fallani, J. E. Lye, V. Guarrera, C. Fort, and M. In- [25 
guscio, Phys. Rev. Lett. 98, 130404 (2007). 

J. Billy, V. Josse, Z. C. Zuo, A. Bernard, B. Hambrecht, [26 

P. Lugan, D. Clement, L. Sanchez-Palencia, P. Bouyer, 

A. Aspect, Nature 453, 891 (2008). [27 

G. Roati, C. D'Errico, L. Fallani, M. Fattori, C. Fort, 

M. Zaccanti, G. Modugno, M. Modugno, M. Inguscio, 

Nature 453, 895 (2008). [28 

C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. [29 
92, 040403 (2004). 

M. Inguscio, W. Ketterle, and C. Salomon (Editors), [30 
Ultra-cold Fermi Gases, Proceedings of the Interna- [31 
tional School of Physics "Enrico Fermi," Course CLXIV, [32 
Varenna 2006, (lOS Press, Amsterdam) 2007. 
R. M. Kalas and D. Blume, Phys. Rev. A 73, 043608 [33 
(2006). [34 

F. M. Cucchietti and E. Timmermans, Phys. Rev. Lett. [35 
96, 210401 (2006). 
K. Sacha and E. Timmermans, Phys. Rev. A 73, 063604 [36 
(2006). [37 
M. Bruderer, W. Bao, and D. Jaksch, Europhys. Lett. 
82, 30004 (2008). [38 
M. Bruderer, A. Klein, S. R. Clark, and D Jaksch, New [39 
J. Phys. 10 033015 (2008). 

D. H. Santamore, and Eddy Timmermans, Phys. Rev. A 
78, 013619 (2008). [40 

D. C. Roberts and S. Rica, Phys. Rev. Lett. 102, 025301 
(2009). [41 
A. Novikov and M. Ovchinnikov, J. Phys. A: Math. [42 
Theor. 42, 135301 (2009). 

J. Tempere, W. Casteels, M. K. Oberthaler, S. Knoop, [43 

E. Timmermans, and J. T. Devreese, Phys. Rev. B 80, [44 
184504 (2009). 

A. Novikov and M. Ovchinnikov, J. Phys. B: At. Mol. [45 
Opt. Phys. 43, 105301 (2010). 

G. D. Mahan, Many-Particle Physics, Plenum Press, [46 
New York 1981. 

F. Chevy, Phys. Rev. A 74, 063628 (2006). 
R. Combescot, A. Recati, C. Lobo, and F. Chevy, Phys. [47 
Rev. Lett. 98, 180402 (2007). 



M. Punk and W. Zwerger, Phys. Rev. Lett. 99, 170404 
(2007). 

N. V. Prokofev, and B. V. Svistunov, Phys. Rev. B 77, 
020408(R) (2008). 

N. V. Prokofev, and B. V. Svistunov, Phys. Rev. B 77, 
125101 (2008). 

M. Punk, P. T. Dumitrescu, and W. Zwerger, Phys. Rev. 
A 80, 053605 (2009). 

A. Schirotzek, C.-H. Wu, A. Sommer, and M. Zwierlein, 
Phys. Rev. Lett. 102, 1 (2009). 

S. Nascimbene, N. Navon, K. J. Jiang, L. Tarruell, M. 

Teichmann, J. McKeever, F. Chevy, and C. Salomon, 

Phys. Rev. Lett. 103, 170402 (2009). 

K. M0lmer, Phys. Rev. Lett. 80, 1804 (1998). 

L. Viverit, C. Pethick, and H. Smith, Phys. Rev. A 61, 

1 (2000). 

S. Yip, Phys. Rev. A 64, 2 (2001). 

R. Roth, Phys. Rev. A 66, 1 (2002). 

H. Pu, W. Zhang, M. Wilkens, and P. Meystre, Phys. 

Rev. Lett. 88, 5 (2002). 

S. Adhikari and L. Salasnich, Phys. Rev. A 76, 1 (2007). 
S. Bhongale and H. Pu, Phys. Rev. A 78, 1 (2008). 
D.-S. Liihmann, K. Bongs, K. Sengstock, and 
D. Pfannkuche, Phys. Rev. Lett. 101, 1 (2008). 
Yu-Li Lee and Yu-wen Lee, arXiv:0910.0603 

B. Ramachandhran, S. G. Bhongale and H. Pu, 

larXiv:0911.2487l 

M. S. Mashayekhi, J. L. Song, F. Zhou.l arXiv:1003.3096l 
T. Karpiuk, M. Brewczyk, S. Ospelkaus-Schwarzer, K. 
Bongs, M. Gajda, and K. Rzazewski, Phys. Rev. Lett. 
93, 100401 (2004). 

T. Karpiuk, M. Brewczyk, and K. Rzazewski, Phys. Rev. 
A 73, 053602 (2006). 

S. K. Adhikari, Phys. Rev. A 72, 053608 (2005). 

G. Bruun, Y. Castin, R. Dum, and K. Burnett, Eur. 

Phys. J. D 7, 433 (1999). 

A. Bulgac and Y. Yu, Phys. Rev. Lett. 88, 042504 (2002). 
M. Grasso and M. Urban, Phys. Rev. A 68, 033610 
(2003). 

A. Niederberger, J. Wehr, M. Lewenstein, and K. Sacha, 
Europhys. Lett. 86 26004 (2009). 

M. Theis, G. Thalhammer, K. Winkler, M. Hellwig, G. 
Ruff, R. Grimm, and J. Hecker Denschlag, Phys. Rev. 
Lett. 93, 123001 (2004). 

Y. S. Kivshar and G. P. Agrawal, Optical Solitons, Aca- 
demic Press, An imprint of Elsevier Science, San Diego, 



California, 2003. 



